On adaptability and "intermediate phase" in randomly connected networks 
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We present a simple model that enables us to analytically characterize a floppy to rigid transition 
and an associated self-adaptive intermediate phase in a random bond network. In this intermediate 
phase, the network adapts itself to lower the stress due to constraints. Our simulations verify this 
picture. We use these insights to identify applications of these ideas in computational problems such 
as vertex cover and K-SAT. 
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The concepts of rigidity and a rigidity transition 0, Q 
have been successfully applied to the study of network 
glasses. The bonds in a network are considered as con- 
straints such that at low connectivity there are more de- 
grees of freedom than constraints, so that the network 
is flexible, or floppy. At high connectivity, there are 
more constraints than degrees of freedom, and the net- 
work is said to be rigid and stressed. The rigidity tran- 
sition lies in between and Thorpe et al. have recently 
suggested the idea of adaptability of a network to avoid 
stress, which has led to the discovery of an intermedi- 
ate phase, between the usual floppy and rigid ones. The 
existence of this phase has been confirmed by numerical 
studies [3 , analysis of finite size clusters 4] and by ex- 
periments 0. What has been lacking is an appropriate 
model of the intermediate phase that allows for analytical 
calculations and provides insight. Furthermore, the in- 
gredients at the origin of the intermediate phase are quite 
generic: adaptability of the underlying network undergo- 
ing a transition, so as to avoid stress. This suggests that 
such intermediate phases may be relevant in many dif- 
ferent fields. For instance, links between rigidity theory 
and computational phase transitions have already been 
suggested 0, but not examined further than an anal- 
ogy. The insight gained from the minimal model we dis- 
cuss here will demonstrate how an intermediate phase 
can arise in computational problems such as the vertex 
cover problem [rf or K-SAT 0. 

We will first present and solve a simple random bond 
model for the rigidity transition. We then introduce 
the possibility of adaptation of the network. Through 
mean-field calculations, we demonstrate the presence of 
an intermediate phase and study its properties. The cal- 
culations are quantitatively confirmed by Monte Carlo 
simulations. We then outline how the results could be 
extended to computational phase transitions. 

Very few non-trivial exactly solvable models are avail- 
able for the standard rigidity transition. To our knowl- 
edge, the randomly bonded models, equivalent to Bethe- 
like lattices 0, are the only ones available. The first 
idea of this paper is therefore to introduce adaptability 
in the randomly bonded models to obtain a description 




FIG. 1: Example of a network with Ni = 11, = 5, so that 



X3 = 0.3125. It has Nn 
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1, Ni3 = 9, Ns-s = 3, iV^ 
N^"' = 3, iV;^"^ = 0, N^"' = 1 (see text). The two subgraphs 
on the right are rigid; the one on the left has three internal 
degrees of freedom: the rotations around the 3 — 3 bonds. 
There is no stress in the network. 



of the complete phase diagram, including the intermedi- 
ate phase. We consider a simple model for the rigidity 
transition. Our network consists of N atoms, A^3 = Nxs 
of which are 3-fold coordinated, and A^i = Nxi are 1-fold 
coordinated; we consider bond stretching as well as bond 
bending constraints. Atoms are bonded randomly, disre- 
garding space: this leads locally to a tree-like network, 
which allows for analytical calculations j^; an example 
with TV = 16 atoms is given Fig. ^ Our aim is not to 
accurately describe network glasses, rather to seek a sim- 
plified description of the salient features. We invoke the 
following constraint counting argument: each atom has a 
priori 3 degrees of freedom (the network is embedded in 
3-dimensional space), each 1-atom brings 1/2 constraint 
(1 bond stretching shared with its neighbor and no bond 
bending), and each 3-atom brings 9/2 constraints (3/2 
bond stretching and 3 bond bending). If no constraint 
is redundant, the number of unconstrained degrees of 
freedom (or fioppy modes) is given by the Maxwell es- 
timate nil] 
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FIG. 2: Comparison of analytical calculations (solid lines) 
with numerics (symbols). The diamonds represent the num- 
ber of floppy modes per atom; the circles (dots) the fraction 
of the network in the percolating rigid (stressed) cluster. The 
dashed line is the constraint counting estimate for the number 
of floppy modes. The numerics are performed on networks of 
4 X 10"' atoms. 



This yields the estimate for the transition point ~ 
0.625, for which the mean number of constraints per atom 
equals 3, the number of degrees of freedom 10] . The cal- 
culation can be carried out exactly using a method de- 
veloped in 11] or by adapting the cavity formalism origi- 
nally devised for spin glasses 12]. The cavity method has 
been used recently 1^ to study a computational phase 
transition: its use is thus the first formal link, beyond 
a simple analogy, between rigidity theory and the field 
of computational phase transitions. These calculations 
confirm the presence of a rigidity transition, first order 
in this case, at a concentration of 3-atoms (i.e. atoms 
with three bonds) — 0.490; see Fig. |21 At the rigid- 
ity transition, a macroscopic rigid cluster appears. In 
the rigid phase, some constraints cannot be fulfilled, and 
stress is present. 

As stress costs energy, it is natural to assume that the 
network will try to adapt itself to avoid it. To make this 
idea more precise, let us define a stress energy for a given 
configuration of the network, equal to the number of re- 
dundant constraints, that is constraints that cannot be 
fulfilled. This energy is of course zero in the floppy phase 
where all constraints can be accommodated, and non- 
zero in the rigid phase, as shown on Fig. 12 Constraint 
counting gives 

Nflop ^3N- iiVi - ^A^3 + Nred , (2) 

where Nfiop is the number of floppy modes in the net- 
work, Nred is the number of redundant constraints, 
3A^ is the a priori number of degrees of freedom, and 
(7Vi-|-9A'3) /2 is the total number of constraints. Thus, for 
a fixed concentration of 3-atoms, counting Nfiop amounts 



to the same as counting N^ed- Consequently, in the fol- 
lowing, we will use indifferently 



H = Nflop otH = N 



red 



(3) 



We assume in addition that the network is adaptive: i.e. 
while Ni and N^ are fixed, some bonds can be rewired to 
decrease the energy thus moving the network away 
from the randomly bonded case. The standard rigidity 
transition described above thus corresponds to the infi- 
nite temperature case. At finite temperature, new phe- 
nomena arise (below). 

To realize the calculations, it is crucial to devise a 
means to measure the degree of organization or random- 
ness of the network. At the crudest level, this is provided 
by the number of bonds between two 1-atoms, one 1-atom 
and one 3-atom, and two 3-atoms, respectively termed 
-^11, A^is, A'33, see Fig.^ In the randomly bonded case, 
a simple calculation shows 
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At finite temperature, Nij can now deviate from its ran- 
dom value N*j, with the constraints 27Vii + N13 = Ni 
and 2iV33 -|- N13 = 3iV3. The level of organization of the 
network, at this crude one bond level, is thus determined 
by a single parameter, which we take to be a = Nn/N^^. 



Fixing a fixes all Ni^ 



1 corresponds to the random 



bonding case, and a 7^ 1 denotes an adaptation of the 
network. Assuming that there are no correlations in the 
network beyond the one bond level, the Chubinsky [lll | or 
the cavity methods allow exact calculations of the energy 
H as a function of X3 and a. It turns out that the lower 
a, the higher is the rigidity threshold for X3. This is in- 
tuitive: decreasing the number of dimers of two 1-atoms 
at fixed X3 decreases the mean connectivity of the main 
network, which is then less likely to be rigid. At fixed 
a, knowing the TVy, it is a simple combinatorial problem 
to compute the configurational entropy of the network; 
S'(x3, a) = Ns{x3, a) = In Q, with 



n(N N N ) - 

il(Nii,Ni3,N33) - 2^1,2^33 Niil N13I N: 
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The conclusion of the calculation now only amounts to 
minimizing with respect to the parameter a the free en- 
ergy H{x3, a) — TS{x3, a), at given T and X3. This step 
is sketched in Fig. |3| What occurs is again intuitively 
clear: for X3 < X3, the random network is unstressed, 
thus entropic and energetic contributions are optimized 
for a — I, and the energy is zero. For X3 > Xg, ran- 
dom bonding has a non-zero energy; however, for small 
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FIG. 3: Energy (dashed line), entropy (dotted line) and free 
energy (solid line) as a function of the parameter a for dif- 
ferent 2:3; from top to bottom, xz — 0.48 (floppy phase), 
X3 = 0.5 (intermediate phase), X3 = 0.53 (rigid phase). The 
circle indicates the rigidity transition varying a at fixed 3:3; 
the diamond indicates the free energy minimum; both coin- 
cide in the intermediate phase. 
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FIG. 4: Comparison between MC simulations (circles for the 
a parameter, diamonds for the number of floppy modes per 
atom) and analytical results (dotted lines correspond to the 
one bond correlation level and solid lines to the two bond cor- 
relation level). The two dot-dashed vertical lines indicate the 
location of the two phase transitions. The dashed line indi- 
cates the Maxwell estimate for the number of floppy modes, 
see eq. Q. MC simulations involve 5000 or 20000 atoms, for 
10*^ MC steps. 



favored. This happens at a second transition X3 = x^* . 
The results are summarized in Fig. 2J To test this pro- 
posed scenario, we performed Monte Carlo simulations 
of the system as follows: at each step, we pick a pair of 
bonds and rewire them, accepting the move according to 
the Metropolis algorithm with energy and tempera- 
ture chosen to be T = 5. Choosing another temperature 
does not induce any qualitative change; only the width of 
the intermediate phase is affected: the lower the temper- 
ature, the wider the intermediate phase. In the limit of 
a non-adaptive network T — > 00, the intermediate phase 
disappears. The network is analyzed, when needed, us- 
ing the pebble game algorithm jlj| , which gives access to 
the number of floppy modes, and thus the energy iPJ, and 
provides the decomposition of the network into rigid and 
stressed clusters. We performed simulations on = 5000 
and N = 20000 systems, for 10*^ Monte Carlo steps. 
Longer test runs did not show any significant differences. 
Results for the parameter a and the number of floppy 
modes, Fig. 01 show a qualitative agreement with the 
simple calculation above: a floppy phase with a = 1 at 
low X3, a stressed phase at high X3 and an intermediate 
unstressed but self-adapted a 7^ 1 phase. 

The agreement is quantitatively not very good: the 
assumption of retaining only one bond correlations most 
likely breaks down. The above calculation can indeed 
be refined to take into account longer range correlations. 
The first steps in this direction consist in counting not 
only the bonds as above, but also the paths along two 
bonds which are of four types: from 1 to 3 to 1, from 1 
to 3 to 3, from 3 to 3 to 1, and from 3 to 3 to 3. This is 
equivalent to counting TVg"-' , N^^^ 



N^^\n^^\ the number 
of 3-atoms that are linked respectively with zero, 1, 2 and 
3 other 3-atoms. Since the two equations 
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must be satisfied, describing the network at the level of 
two bond correlations requires the introduction of two 
new parameters, defined as follows: 
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enough 2:3, decreasing a to avoid stress energy costs lit- 
tle entropy, and is favorable. This is the intermediate 
phase. Throughout this phase, the system adapts itself 
to stay exactly on the verge of stress, so that the en- 
tropic cost is minimal and there is no stress. For X3 
large enough, decreasing a further to avoid stress costs 
too much entropy, and the appearance of some stress is 



where N^^"^* is the number iVg*'' expected without two 
bond correlations. Thus, any 7^ 1 shows the pres- 
ence of two bond correlations. Using the cavity method, 
one can calculate the energy 7?(x3, a, ao, ai). As above 
with the parameter a, it turns out that the lower ao and 
ai are, the less stress-prone is the network. Evaluating 
the entropy S{x3,a,ao,ai) is still a simple combinato- 
rial problem. The results are plotted in Fig. ^ and are 
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FIG. 5: Diamonds (circles): probability that a rigid (stressed) 
cluster percolates the entire sample, as a function of x^. The 
runs were made with N = 5000 or A'' = 20000 atoms, for 10*^ 
MC steps. Solid and dashed lines are guides for the eyes; dot- 
dashed lines indicate the location of the two phase transitions. 

much closer to the numerical results. Due to the simplic- 
ity of the model, we were able to make the calculations 
up to three bonds correlations; they confirm the trend 
toward the numerical results, but do not completely co- 
incide with them in the intermediate phase. This leads 
to the important conclusion that medium or long range 
correlations are important in the network, especially in 
the intermediate phase j^] . 

Monte Carlo simulations also provide the opportunity 
to examine the different phases in greater detail. The 
most important question concerns the presence or ab- 
sence of percolating rigid or stressed clusters. In the 
floppy phase, there is no percolating rigid or stressed 
cluster, whereas in the stressed phase rigidity and stress 
always percolate. The interesting case is the intermedi- 
ate phase: our simulations show that the probability to 
find a rigid or stressed percolating cluster varies smoothly 
from to 1, see Fig.|31 The probability to find a stressed 
cluster is always slightly lower than the probability to 
find a rigid one. Thus, while in the intermediate phase 
the system lies right on the boundary between rigid and 
floppy, the probability to be on one side or another of the 
boundary goes from to 1. Although strictly speaking 
the model gives access only to Monte Carlo sampling of 
the phase space, we can attempt to infer a dynamical pic- 
ture of the intermediate phase, namely a percolating rigid 
cluster disconnecting and reconnecting elsewhere through 
local connections. 

From the detailed analysis above, the ingredients 
leading to the intermediate phase appear to be clear: 
an underlying first order phase transition, and the 
possibility of adaptation of the network, giving rise 
to an entropy competing with the energy associated 
with the underlying transition. One thus expects the 



same phenomena for a variety of computational phase 
transitions. Consider the vertex cover problem as an 
example. The problem is as follows: given a network 
of size N, is it possible to find a subset of vertices of 
size xN, such that all edges are adjacent to this subset? 
Such a subset is called a vertex cover. If the network 
structure is random and fixed, this is always possible 
at low connectivity, and becomes impossible at high 
enough connectivity 0, through a phase transition. If 
one defines the energy of a configuration as the number 
of edges not adjacent to the subset chosen, and allows 
for a reorganization of the network, an intermediate 
phase with a finite probability of finding a vertex cover is 
likely to appear. This should apply also to satisfiability 
problems of the K-SAT type . More generally, as most 
networks found in physics, social science or biology, are 
adaptive, the intermediate phase described in this work 
is likely to be encountered in a wide variety of situations. 
Future work will elaborate the thermodynamic and dy- 
namical signatures of the intermediate phase, including 
expectations of filamentary geometry and glassiness. 



We would like to thank Don Jacobs, Mykyta Chu- 
binsky and Michael Thorpe for permission to use the 
software Pebble-3D. Motivational discussions with Jim 
Phillips are gratefully acknowledged. Work at Los 
Alamos National Laboratory is supported by the US De- 
partment of Energy. 



[1] J. C. Phillips, J. Non-Cryst. Solids 34, 153 (1979). 
[2] M. F. Thorpe, J. Non-Cryst. Solids 57, 355 (1983). 
[3] M. F. Thorpe, D. J. Jacobs, M. V. Chubinsky, and 

J. C. Phillips, J. Non-Cryst. Solids 266-269, 859 (2000). 
[4] M. Micoulaut and J. C PhiUips, Phys. Rev. B 67, 104204 

(2003). 

[5] D. Selvenathan, W. J. Bresser, and P. Boolchand, Phys. 

Rev. B 61,15061 (2000). 
[6] R. Monasson, R. Zecchina, S. Kirkpatrick, B. Selman, 

and L. Troyansky, Nature 400, 133 (1999). 
[7] M. Weigt and A. K. Hartmann, Phys. Rev. Lett. 84, 6118 

(2000). 

[8] M. Mezard and R. Zecchina, Phys. Rev. E 66, 056126 
(2002). 

[9] P. M. Duxbury, D. J. Jacobs, M. F. Thorpe, and 
C. Moukarzel, Phys. Rev E 59, 2084 (1999). 

[10] This argument can actually be refined as shown in [Till : 
it then gives a much better result. 

[11] M. V. Chubinsky, PhD Thesis (Michigan State Univer- 
sity, 2003). 

[12] J. Barre, A. R. Bishop, T. Lookman, A. Saxena, in prepa- 
ration. 

[13] M. Mezard and G. Parisi, J. Stat. Phys 111, 1 (2003). 
[14] D. J. Jacobs and M. F. Thorpe, Phys. Rev. Lett. 75, 
4051 (1995). 



